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QQ Abstract 

I— I The quantum entanglement dynamics of a one-dimensional spin-1/2 anisotropic XXZ model is 

studied using the method of the adaptive time-dependent density-matrix renormalization-group 
C when two cases of quenches are performed in the system. An anisotropic interaction quench and 
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the maximum number of domain walls of staggered magnetization quench are considered. The 
dynamics of pairwise entanglement between the nearest two qubits in the spin chain is investigated. 



I> The entanglement of the two-spin qubits can be created and oscillates in both cases of the quench. 

The anisotropic interaction has a strong influence on the oscillation frequency of entanglement. 

^ PACS numbers: 03.67.Mn, 03.65.Ud, 75.10.Pq 
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I. INTRODUCTION 



Entanglement generation and distribution is an important problem in performing quantum- 
information tasks, such as quantum computation and quantum teleportation [TH3]. Many 
results showed that entanglement existed naturally in the spin chain when the temperature 
is at zero. Due to its potential applications, the pairwise entanglement of anisotropic Heisen- 
berg model was extensively studied recently [3]. A typical example is the exactly solved 
one-dimensional infinite-lattice anisotropic XY model using the Jordan- Wigner transform 
to obtain the pairwise entanglement [5H8]. By using the Bethe ansatz solution, the entan- 
glement in infinite isotropic Heisenberg rings was calculated fO]. In recent years, the study 
of the dynamics of entanglement has attracted much attention due to the manipulation of 
quantum systems. The dynamics of entanglement in different systems were investigated 
|10t4T9] . The effects of anisotropy interaction on the evolution of entanglement were investi- 
gated [I9] , such a global quench could actually be realized for atoms in optical lattices [20] . 
Moreover, the effects of a bond defect on the entanglement dynamics in a local quench were 
studied [2T] and the evolution of one-dimensional quantum lattice systems was investigated 
when the ground state was perturbed by altering one site in the middle of the chain [221 - I27] . 
Recently, Bose exploited a global quench dynamics in spin chains for distant pairwise en- 
tanglement, which could be used for quantum communication [28]. It would be interesting 
to investigate the overall pairwise entanglement dynamics between the nearest two qubits 
in the spin chain when several typical quenches are performed in the system. 

In this paper, quench dynamics of entanglement in an opened anisotropic Heisenberg spin 
chain is analyzed. In section II, the Hamiltonian of an opened anisotropic Heisenberg spin 
chain is presented. In section HI, the dynamics of pairwise entanglement between the nearest 
two qubits in the spin chain is studied when two kinds of typical quench are performed in 
the system. In section IV, a discussion concludes the paper. 

II. HAMILTONIAN OF THE SYSTEM 

The Hamiltonian of an opened Heisenberg XXZ chain with N sites is given by 

N-l 

H = $^[Jx.,(ara|;, + afaf+J + J.a,V,^J, (1) 

i=l 
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where erf (a = x,?/, z) are Pauli operators on the i-th site, is the length of the spin chain, 
Jxy, Jz denotes the couphngs in the a;?/-plane and the 2;- axis respectively. For simplicity, 
the couplings in the plane J^y = 1 is considered. An opened boundary condition (OBC) is 
assumed because the antiferromagnetic Heisenberg spin chain with OBC can be achieved 
artificially in the experiment |29j . 

In the paper, the concurrence is chosen as a measurement of the pairwise entanglement 
[3]. The concurrence C is defined as 

C = max{Ai - A2 - A3 - A4, 0}, (2) 

where the quantities Aj(i = 1, 2, 3, 4) are the square roots of the eigenvalues of the operator 
Q = Pi2(cri ® cr|)pi2(o"i o"f ). They are in descending order. The case of C = 1 corresponds 
to the maximum entanglement between the two qubits, while C = means that there is no 
entanglement between the two qubits. 

III. PAIRWISE ENTANGLEMENT DYNAMICS 

In this section, the entanglement dynamics of a spin— 1/2 antiferromagnetic Heisenberg 
chain is analyzed when two kinds of quenches are performed in the system. One kind of 
quench is performed by abruptly varying the anisotropic interaction from very large value 
down to finite values. Another kind of quench is performed by sudden release of the domain 
walls of staggered magnetization. It is known that it is hard to calculate the dynamics 
of entanglement because of the lack of knowledge of eigenvalues and eigenvectors of the 
Hamiltonian. For models that are not exactly solvable, most of researchers resort to exact 
diagonalization to obtain the ground state for small system size. While for small system size, 
some phenomena might be missing. If the system size is large, other kinds of phenomena 
can appear and cen be seen much clearer. 

For large system size, the adaptive time-dependent density-matrix renormalization-group 
(t-DMRG) can be applied with a second order Trotter expansion of Hamiltonian as described 
in [30l El]. In order to check the accuracy of the results of t-DMRG, the results of exact 
diagonalization can be considered as a benchmark for a small size system. In the simulation, 
a Trotter slicing 6t = 2.5 ■ 10~^ and Matlab codes of t-DMRG with double precision are 
performed with a truncated Hilbert space of m = 60. It turns out that a typically discarded 
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weight of Sp = 10 ^ can keep the relative error SC in C below 10 ^ for a chain of L = 60 
sites with time t < 40/ J. 

A. Anisotropic Interaction Quench 

To start t-DMRG, some initial states need to be prepared. Firstly, the initial state can be 
prepared to be a perfect antiferromagnetic Neel state. The perfect antiferromagnetic Neel 
state is the ground state of Eq. (1) with — ?■ +oo. It is given by 

|zm)i = |titi---n)- (3) 

Such a state has been achieved with high fidelity by using decoupled double wells |32]- 
The procedure adopted for the initial state \ini)i is to calculate it as the ground state of 
a suitably chosen Hamiltonian Hint = H{Jz — +oo). The quantum quench is performed 
when the anisotropic interaction parameter is abruptly varied from an initially very large 
value down to finite values [28] . 

The entanglement Cj^j+i between the nearest two qubits is plotted in Fig. 1 as a function 
of the spin site i and the time t when = 1.0. The three dimensional entanglement Cj_j+i 
is plot in Fig. 1(a). The contour line is plotted in Fig. 1(b). From Fig. 1, it is clear 
that there is an entanglement wall for small time t. It seems that the sudden change of 
Jz from extremely large value (J^ — )■ +oo) to a finite value creates the entanglement wall 
[71 IE]. The entanglement in the range of 10 < i < 50 is quite similar. When the time t = 0, 
the entanglement is zero. When the time increases, the entanglement creates and increases 
from zero. The entanglement reaches a peak at the time t = 0.7. When the time increases 
further, the entanglement value maintains in the extent [0,0.1]. In Fig. 1, it could also be 
seen that when i < 10 and 50 > i, the entanglement oscillates acutely. The maximal value 
of entanglement is higher than the maximal value in the extent of 10 < i < 50. This is 
due to the effects of the open boundary condition [191 EH]- From Fig. 1, it is clear that the 
entanglement Cj_j+i is symmetric about i = 30. This is mainly due to the fact that several 
sources in a lattice emit oppositely moving pairs of entangled quasi-particles [8]. 

In order to avoid the boundary effect, the entanglement Cs^si of the central two qubits 
is plotted as a function of the time t for different anisotropic interaction in Fig. 2. The 
anisotropic interaction has little influence on the height of the first peak of the entangle- 
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ment. When = 0, 0.5, the entanglement drops down to zero after reaching the peak. 
The entanglement cannot generate at any longer time. When = 1.0, the entanglement 
can appear and disappear alternatively. When = 1.5,2.0, the entanglement can create, 
oscillate and does not disappear for longer time. Except the first peak in Cs^si, the height 
of oscillation of entanglement with Jz = 2.0 is higher than that with Jz = 1.0, 1.5. This 
is due to the strong dependence of the damping coefficient of the entanglement wave on 
the anisotropic parameter. When t increases, the relative height of the first peak and the 
average value of C3o,3i becomes smaller with increasing J^,. The increase of Jz reduces the 
relative height in 6*30,31 [11 [S] . 

B. Domain Walls of Staggered Magnetization Quench 

Secondly, the initial state can be prepared to be an inhomogeneous initial state. It can 
be given by [3H |35] 

\ini)2 = I ti • ■ ■ tN/2iN/2+i ■■■In), (4) 

where all spins on the left half are pointing up along the z axis, while all spins on the right 
half pointing down. The state contains many high-energy excitations and is thus far from 
equilibrium. It can be considered as a state with almost the maximum number of domain 
walls of staggered magnetization. The quantum quench is performed when the domain walls 
are suddenly released in a Heisenberg XXZ chain with different anisotropic interaction [35] . 

The entanglement Cj,j+i between the nearest two qubits is plotted as a function of spin 
site i and time t in Fig. 3 when Jz = 0. The three dimensional plot of Cj,i+i is shown in 
Fig. 3(a). The counter line of Cj^j+i is plotted in Fig. 3(b). From Fig. 3, it is seen that 
the entanglement Ci.i+i can be generated. The entanglement is mirror symmetric about the 
line of 2 = 30. The entanglement 6*30,31 generates firstly, then the entanglement 6*29,30 and 
6*31,32 generate, and so on. The period of entanglement 6*j,j+i generation is a linear function 
of |i — 30|. The notion of diffusive dynamics of pairwise entanglement in the system can be 
clearly seen [351 EE]. The two central neighboring sites are in the state | tN/2iN/2+i)- Due 
to the exchange interaction, they will be entangled after the time t = 0. The entanglement, 
initially localized on the two central neighboring sites of the chain will spread owing to the 
exchange interaction [TJ |8] . 
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In order to see the effects of the anisotropic interaction, the entanglement 6*30,31 and 
the entanglement 6*25,26 of the central two qubits is plotted in Fig. 4 as a function of the 
time t for different anisotropic interaction. In Fig. 4(a), the entanglement 6*30,31 generates 
immediately after the time t = 0. The entanglement reaches the first peak rapidly and then 
oscillates around a constant. The anisotropic interaction has a weak influence on the first 
peak of the entanglement. The first peak decreases slightly when the anisotropic interaction 
increases. When = 0.0, the entanglement oscillates around 0.35. When = 0.5, the 
entanglement oscillates around 0.12. When = 1.5, 2.0, the entanglement oscillates around 
0.40 and exists all the time. For the above cases, the entanglement maintains for a long 
time. While for = 1, the entanglement appears and disappears alternatively. After 
t = 37, the entanglement disappears all the time. It is clear that the entanglement 6*30,31 of 
the central two qubits decreases as the anisotropic interaction Jz increases and then increases 
as Jz increases further and finally almost saturates. In Fig. 4(b), the entanglement 6*25,26 
generates after a short time. This time delay is caused by the finite propagation speed of 
the entanglement from the central pair to side pairs. The entanglement reaches the first 
peak rapidly and then oscillates around a constant. The anisotropic interaction has a strong 
influence on the first peak of the entanglement. The height of the first peak decreases 
when the anisotropic interaction increases. When Jz = 0, the first peak is 6*^5 26 = 0.427, 
while the peak is 6*25 26^^^ = 0.259 when Jz = 0.5. The anisotropic interaction has a feeble 
influence on the period when the entanglement reaches the first peak. The time of the 
entanglement reaching the first peak decreases slightly when the anisotropic interaction 
increases. After the entanglement reaches the first peak, the entanglement oscillates around 
0.33 when Jz = 0.0. When Jz = 0.5, the entanglement oscillates around 0.11. For the 
above cases, the entanglement exists all the time. While for Jz = 1, the entanglement 
appears and disappears alternatively. After t = 25, the entanglement disappears all the 
time. When Jz = 1.5,2.0, the entanglement disappears and never appears after the first 
peak emergences. The entanglement decreases as the anisotropic interaction increases. It 
is noted that the propagation velocity of the entanglement is slightly influenced by the 
anisotropy parameter, but the damping coefficient of the entanglement wave is strongly 
dependent on the anisotropic parameter^ |8]. . 
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IV. DISCUSSION 



By using the method of the adaptive time-dependent density-matrix renormahzation- 
group, the time evolution of the entanglement in a one-dimensional spin-1/2 anisotropic 
XXZ model is investigated when two quenches are performed. One quench is a sudden 
change of the anisotropic interaction, while another is the abrupt release of the maximum 
number of domain walls of staggered magnetization. The dynamics of pairwise entanglement 
between the two nearest qubits in the spin chain is studied. The entanglement of the two- 
spin qubits can be created and oscillates after one of the two quenches is performed. In both 
cases of the quench, the time evolution of the entanglement is symmetric about the central 
qubit. For the anisotropic interaction quench, an entanglement wall creates. Except the first 
peak, the entanglement increases and oscillates as the anisotropic interaction increases. 
For small J^, the entanglement may disappear at longer time. In the quench of domain 
walls of staggered magnetization, the entanglement of central two qubits appears first, after 
a short period, the entanglement of the next two qubits appears, and so on. The height 
of the first peak in the entanglement also decreases as the pair of two qubits moves away 
from the center. The entanglement of the central two qubits decreases when the anisotropic 
interaction increases, and then increases when the anisotropic interaction increases further. 
For two kinds of the quenches, the oscillation frequency of the entanglement is strongly 
infiuenced by the anisotropy. This phenomenon may be used to control the dynamics of the 
entanglement by varying the anisotropic interaction of the Heisenberg spin chain. 
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FIG. 1. The pairwise entanglement Cj^j+i between the two nearest qubits is plotted as a function 
of spin site i and time t with = 1 when there is anisotropic interaction quench. The size of the 
system is = 60. (a). The three dimensional plot of Ci^i-^-i. (b). The contour line of Cj^j+i. 
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FIG. 2. The pairwise entanglement Cao.ai of the two central qubits is plotted as a function of time 
t for different anisotropic interaction when there is anisotropic interaction quench. The sizes of 
the system is = 60. 
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FIG. 3. The pairwise entanglement Cj^j+i between the two nearest qubits is plotted as a function 
of spin site i and time t with = 1.0 when there is domain walls quench. The size of the system 
is = 60. (a). Three dimensional plot of (7j_j+i. (b). The contour plot of (7j_j+i. 
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FIG. 4. The pairwise entanglement is plotted as a function of time t for different anisotropic 
interaction Jz when there is domain walls quench. The size of the system is = 60. (a). The 
pairwise entanglement C25,26- (b). The pairwise entanglement 6*30,31 of the two central qubits. 
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